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I. INTRODUCTION 



The factorization method for solving the eigenvalue equations for the Schrodinger operators was, first introduced 
by Schrodinger which he also applied it to the hypergeometric differential equation Later, the method was 
generalized by Infeld and Hull , and an extensive list of references of the method can be found in Q ■ It is an effective 
algebraic method which is related to the mathematical structure of supersymmetric quantum mechanics 0-0| and 
also to the concept of the shape invariance . 

The linear harmonic oscillator is the textbook-example for which the method is applied elegantly. In this case, all the 
analytical properties of the eigenfunctions, namely, the Hermite polynomials, can be derived from the factorization 
method itself. However, this is not the case for the multi-dimensional isotropic harmonic oscillator for which the 
special function related to the radial eigenfunctions is the generalized Laguerre polynomials. The generalized Laguerre 
polynomials occur in various other well-known quantum mechanical problems, such as the Morse oscillator and the 
Coulomb problem and thus their properties are essential for the study of these problems as well. 

One of the algebraic methods to study the spectrum of the isotropic harmonic oscillator problem is to employ a 
realization of so (2, 1) Lie algebra as a spectrum generating algebra [13] ■ Another algebraic method is to transform 
the radial equation into a confluent hypergeometric equation and then factorize the resulting equation which requires 
the use of the properties of the generalized Laguerre polynomials [3, 4]. Still another algebraic way is to start from 
the generalized Rodrigues formula for the Laguerre polynomials and define ladder operators from various recurrence 
relations [TTI . [T2j which again requires the use of the properties of the Laguerre polynomials. 

In the present paper, we will show that the properties of the generalized Laguerre polynomials can be obtained as 
a by product of the factorization of the radial equation for the isotropic harmonic oscillator in arbitrary dimensions. 
We also show that the same procedure can be applied to obtain bound energy eigenstates of the Morse oscillator as 
well. 

The paper is organized as follows. In Section II we set our notation by defining the ladder operators that factorize 
the radial part of the Schrodinger equation for the isotropic harmonic oscillator, following the method of Infeld-Hull 
Using these ladder operators, we also define new type of Z-changing ladder operators which leave value of radial 
quantum number j fixed. We also show how to construct further new ladder operators that shift between any given 
pair of radial eigenstates labeled by angular and radial quantum numbers 

In Section HI we indicate how the properties of the generalized Laguerre polynomials can be obtained by using the 
various ladder operators and the radial eigenfunctions defined in Section II. 

In Section IV we find the energy eigenstates of the Morse oscillator in terms of the generalized Laguerre polynomials 
by using the ladder operators obtained by factorizing the Morse hamiltonian. In Section V we present our conclusions. 
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II. ISOTROPIC HARMONIC OSCILLATOR IN d-DIMENSIONS 



The details of the factorization method presented here can be found in Q where the isotropic harmonic osciUator is 
of classified as type-C factorization. We also refer to 01 for the factorization method applied to the radial Schrodinger 
equation for a non-relativistic particle with the central potential of the form ^muj^r^ + in d dimensions from a 

different point of view. 

We shall work with the variables expressed in natural units. The natural unit of energy for the harmonic oscillator 
is huj. Dividing the Schrodinger equation by huj and also scaling length by rg = momentum by H/vq one 

has the Schrodinger equation in terms of the natural units. In the natural unit system, the canonical commutation 
relations read [Xm,Pn] = i5mn where m,n = 1, 2, ... d. 

Up to a sign, square of the momentum in coordinate representation is given by the generalized Laplacian operator. 
In Euclidian space K'', and in terms of spherical coordinates, it can be written as 

Am"* 1 d J T d 1. qd-l ,^ ^ 

A = r \ A (1) 

where A'^ is the Laplacian operator defined on the unit sphere S'^~^ with the metric induced from that of the 
ambient Euclidian space IR'^. The radial momentum, in accordance with the Weyl prescription of ordering of non- 
commuting operators, can be defined as 

in the position representation. The radial momentum has the canonical commutator [r,^^] — i with the radial 
coordinate and one has Pr — pi with respect to the inner product defined in terms of the the weight function 
w{r) = r^'^~^^ with r G [0,cxi). With this weight function, the inner product of two functions |/) and |g), in the 
Hilbert space of functions spanned by the square integrable energy eigenfunctions becomes 

/•oo 

{f\g} = / drr''-'r{r)9{r). (3) 



r- -_=p^ + (d_l)(d-3)^. (4) 



Jo 

The radial term of the Laplacian operator in ([l} can be written in terms of the coordinate representation of the radial 
momentum as 

(d — 1) number of angular coordinates can be separated from the radial coordinate by introducing generalized spherical 
harmonics, yi{f), in d dimensions fl3|. Recalling that r'3^;(f) can be expressed in terms of homogeneous irreducible 
monomials of order I in the Cartesian coordinates Xm- Hence, they satisfy A^ r'^yiif) = and therefore using ^ 
one has 

^^"~"yi{f) = -i{i + d-2)yi{f). (5) 

As a result, one obtains the set of effective hamiltonians hi, depending on the radial coordinate and also on the 
separation parameters I and d. hi can be written as 



2 
Pr 



This form of hi implies that the eigenvalue equation hi\l,i) — ej{l)\l,j) admits factorization since it is of the form the 
sum of the squares of two operators up to a constant remainder. Note also that the radial Coulomb-type potential of 
the form ^ also allows the corresponding effective radial hamiltonian to be written as the of product of two hermitian 
conjugate operators up to a constant remainder. 

The first label of radial eigenkets \l,j) refers to the angular momentum quantum number, whereas the second 
one j refers to the radial quantum number. For convenience, the energy eigenvalues corresponding to the radial 
eigenfunctions i?Q''(r) —: {r\l,0) are denoted by ej{l), all of which depend on the parameters l,j and d. 

The ladder operators, which are the conjugates of each other with respect to the inner product ([3]) defined with 
respect to the weight function w{r) = r'''"^-', are of the form 

Di = ^(+zpr-{2l + d-l)^+ry (7) 
Dj = ^(^~ipr-{2l + d-l)^+ry (8) 
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By using the same symbols for the position representations for the ladder operators, the position representations of 
these operators can be written as 

A - +^^'e-^V.|.,.VV-, (9) 

The particular forms of the position representations ([9]) and ([TOl) for the first order ladder operators can be obtained 
from ([7]) and ^ respectively by introducing appropriate integrating factors for the multiplicative terms of the ladder 
operators and will be shown to be very useful for the discussions below. 

For I I—)- and d i— >■ 1, the expressions in ^ and pO)) become the position representations for the ladder operators 
of the one dimensional harmonic oscillator, namely a and respectively. 

The choice of the ladder operators having the same algebraic structure is not unique. However, the eligible ladder 
operators can be obtained from the requirements: (i) The key functions which are annihilated by Di should be 
normalizable, that is, {1, 0\l, 0) = 1, (ii) Sj{l) < ej+i{l) for all /, j = 0, 1, 2, . . . which ensures the normalization of the 
radial eigenkets that is, = 1 for j > 0. In terms of the ladder operators above, the hamiltonian can be 

written as 

hi=D\Di+e{l) (11) 

where eo(0 = I + d/2 is the ground state energy eigenvalue corresponding to the key function i?Q provided that 
DiB!^^ — 0. The spectra of the sequence of the hamiltonians h(^i^j) with /, j = 0, 1, 2 . . . follow from the algebra of the 
ladder operators (see Fig. 1). From the definitions of the ladder operators, one finds the recurrence relations 

+ 2 - (12) 
with Z, J = 0, 1, 2 . . .. In terms of the sequence of the hamiltonians which are naturally defined as 

^C+j) = + £j(0- (13) 

The recurrence relation p2p for j = can be written as 

hlDj=Djh^l+,y (14) 

Iteration of this result in the radial quantum number j gives 

The recurrence relation implies that the ground state energy of h^i^j) corresponds to the j*'* excited level of 
hi which is given by ej{l) — I + 2j + d/2. Thus, the radial energy eigenfunction corresponding to the j"* excited 
eigenstate of the hamiltonian is given by 

i^'j) «^i^i+i)^i+2)---4+.-i)i^+^'0) (16) 

up to an appropriate normalization constant to be found below. The radial key eigenfunctions are annihilated by the 
ladder operator 

i?(,+j)|Z + j,0) =0. (17) 
This yields the normalized key radial eigenfunctions as 

{r\l + J, 0) = 4'+^' (r) = [ir(Z + J + d/2)] "^/^ r^'+^^e-^'/' (18) 

by inspecting the position representation of the ladder operator in ([5]). It follows from both the definition of the key 
functions and of the ladder operators in harmony that 

i?(,+j)|Z + j,0) =ri?(,+,_i)r-i|Z + j,0) =0 (19) 



hiDjD^,^^, ■ ■ ■ d].^^_^, = d\d].^^, ■ ■ ■ D],^^_^Mn+^) . (15) 
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and as a result one finds |/ + j, 0) oc r^\l + j T IjO). Therefore, the operators r^^ can be identified as the ladder 
operators which shift between the successive key functions. Furthermore, the definition of the operator in (|10P 

allows the ladder operators to be written in terms of Dq as 



(20) 



The relation in and its conjugate are very practical in the calculations below. Using ^U\\ . it is possible to rewrite 
all the excited radial eigenstates of the hamiltonians hi as 



\l,j) (X DlDl^,^Dl^^yDl^^_,^\l + j,0), 

1 d^' 



^[2(/+i-l)+d]g-r^ 



(21) 

(22) 



up to a normalization constant. The normalization of the eigenfunctions Rj (r) follows from the normalization of the 
eigenfunctions i?Q "''"'^ (r). Using and by induction one finds 



n(£fc(0-eo(0)('+J,0|?+j,0) 



fc=i 



= 2^j!(/+j,0|/+j,0). 
Therefore, the normalized radial eigenfunctions can be written in a convenient form as 



Rf\r) = 2(^'-i)j!r(/+j+d/2) 



-1/2 



d(r2) 



^2(;+rf/2-i+i)g- 



(24) 
(25) 



(26) 



The expression p6p is a generalized Rodrigues-type formula for the radial eigenfunctions and it is easy to write it in 
terms of the generalized Laguerre polynomials in the variable as 



^ ' ^ ' [T{l+J+d/2)J '^'^^ ' 



-{l+d/2-l) 



(27) 



{z) defined 



for J > 1 values of the radial quantum number. In the next section we will show that the functions 
above do indeed satisfy the generalized Laguerre differential equation in the variable z = by using the recurrence 
relations provided by the ladder operators, cf. Eqns. ([57)) - (H(Il) below. Note that the radial key eigenfunctions 
themselves are not generalized Laguerre polynomials and that the minimum value of I + d/2 — 1 is —1/2 which 
corresponds to values of ^ = and d = 1. The Laguerre polynomials, namely, 's, occur only for d = 2 and zero 
angular momentum radial energy eigenstates. The expressions ((26)) and ((27)) hint at the possibility that the ladder 
operators can be used to obtain the properties of the generalized Laguerre polynomials which are usually obtained by 
other means [m.The functions e-''/^zf'L'^{z) are sometimes called the generalized Laguerre functions [12| . In Section 
IV, we will also show that for z = 2e~^ these functions correspond to the bound eigenstates of the Morse oscillator. 
In the study of the coherent states defined in terms of the ladder operators, a hamiltonian whose eigenfunctions 
are generalized Laguerre functions were constructed in (T^). In a more general scheme, it is also possible to define 
ladder operators for the classical polynomials starting from the corresponding generalized Rodrigues formula and the 
recurrence relations [H, [l^ . 

It is possible to construct a new type of ladder operators, namely ^-changing ladder operators, which shift between 
the radial eigenfunctions with a fixed value of the radial quantum number j. For this purpose, first recall that 
\l + j ^ 1,0) oc r^^\l + j, 0) and therefore, using this property, the relation (IT6)) can be written as 



D 



(;+i)-^(;+2)-^(;+3) 



'D]D 



■Dl^,yr^\l+m 



(28) 



up to a normalization constant. By commuting to the left of the product of the ladder operators, it is possible to 
rewrite the right hand side back in terms of \l,j) again. To this end, one needs the commutator identity 



D 



■ r^ 



= -V2jrD 



(i+l)^(/+2) 



■D 



{l+]-2)^(l+j-l) 



(29) 



£o(0 + 6 



eo(0 + 4 



eo(0 + 2 



soil) 



Dj |Z + 1,2). Am) 
\l,2) -Cl(;^i)|/ + 2,1) ^D, 



/ + 2,2) : K + 4,1) 

I' + 3,1) 



|/ + 5,0> 



7 ic^ J^(l+2) 



|; + 4,o> 



(i+2)|/ + 3,0) 



|Z,1) -^^'-^^ |/ + 2,0) 

\l,0) 



I 



l + l 1 + 2 1 + 3 1 + 4: 



1 + 5 



FIG. 1: The lattice of the energy levels for the sequence of hamiltonians h{i^jy The column of the equally spaced eigenkets 
above the ground level \l + j, 0) for each j = 0, 1, 2 . . . constitute the spectrum of the hamiltonian h^i^jy The energy difference 
between successive key eigenkets is one unit (of &j) whereas the energy difference of the two successive eigenkets of the 
hamiltonian is two units. The actions of the ladder operators on the eigenkets are indicated by the arrows. 



The commutator (|29p can be found by using (j20p and the canonical commutation relation [r, p,.] = i. Using this result 
(Pg)) takes the form 

|/ + (X -V23\l + l,i - 1) + r|/,j) (30) 
up to an overall constant. By using (|26p . this can be written in terms of the ladder operators as 

|/ + l,j)(xe'-'Ae-'-'|/,j). (31) 

again up to an overall constant. The operator on the right hand side of ([5l]) evidently defines an Z-changing ladder 
operator which will be denoted by . Using the above notation, the action of this new ladder operator on the 
radial eigenkets can be written as 

Vljl,j)^\l + l,3). (32) 

The conjugate of this operator defines yet another Z-changing ladder operator, j), which can be simply written as 

\l + 1, j) (X I?(zj)|Z, j) = e'-"Dy\l,j) (33) 

up to a constant depending on both of the quantum numbers. These new Z-changing ladder operators ^-^ and T^{i,j) 
depend on both of the quantum numbers I and j and also on the dimension d which is not indicated for convenience 
of the notation. Note here that, it is possible to define an operator composed of appropriate combinations of the 
ladder operators defined so far, namely, -^Ci) ^'^'^ t"^^, that shifts between any given pair of the 

eigenkets and iUustrated in Fig. 1. For instance, the combination DjvJ^^ .-^ of the ladder operators acting 

on the eigenstate \l,j) gives 

Djvl,Jl,j)^\l,j + l) (34) 

up to a normalization constant. Therefore, it is possible to identify ojiy^ and its conjugate, namely D^i ffDi, as 
j-changing ladder operators. 

It is worth to emphasize at this point that (i)The ladder operators above can be used to construct new ladder 
operators that shift between any given states \l,j) and \l',j') and Z-changing and j-changing operators defined above 
are examples of these new of ladder operators which shift between |Z + l,j) and \l,j + 1) respectively, (ii) Nowhere in 
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the above construction the properties of the generahzed Laguerre polynomials have been used. Furthermore, various 
recurrence relations for the generalized Laguerre polynomials can be derived using any of the ladder operators defined 
above. From this point of view, the presentation above can be considered to unify the method of solving eigenvalue 
problem and obtaining the recurrence relations of the related special function as a by product. In contrast to the 
other approaches, the ladder operators are constructed by invoking the properties of related special functions 
and the matrix elements for any of the ladder operators defined above can be calculated using the mathematical tools 
constructed above without employing any of the properties of the Laguerre polynomials as additional mathematical 
property of the solutions. For example, integrals involving the generalized Laguerre polynomials can be calculated by 
using the normalization of the radial eigenfunctions and the expressions of the radial eigenfunctions in terms of the 
ladder operators, in particular, the conventional normalization of the generalized Laguerre polynomials can be found 
from the normalization of the radial eigenfunctions. 

As will be shown in the next section, the expressions involving the ladder operators and the energy eigenfunctions 
can be translated into the expressions involving the generalized Laguerre polynomials. Since it is possible to construct 
a number of different recurrence relations for the radial eigenfunctions by defining an appropriate composition of the 
ladder operators defined above, the factorization method also provides a new way of deriving new recurrence relations 
for the generalized Laguerre polynomials as well. 



III. PROPERTIES OF THE GENERALIZED LAGUERRE POLYNOMIALS 

For the definitions and the conventions of the Laguerre polynomials adopted here we refer to Pill. First, we note 
that the standard normalization of the generalized Laguerre polynomials follows easily from the normalization of the 
radial eigenfunctions, all of which can be turned into Gamma function-type integrals. The radial eigenfunctions which 
are normalized to unity are given explicitly in (j26p and they can be expressed in terms of the Laguerre polynomials 
using (1271) . Therefore, the normalization of the radial eigenfunctions \l,j) can be expressed in terms of the Laguerre 
polynomials as 

/•OO 

(0/ 



nj + 1) 



T{l + d/2 + j) 

where the variable = z has been introduced in the second line. From this expression it follows that 

'dze-zn£,^(z)|^ = '^^^+j+'^ (35) 

with ^ — I -\- d/2 — 1. This standard normalization can also be obtained by means of the Rodrigues formula for the 
Laguerre polynomials. 

Second, we illustrate that any expression for the radial eigenfunctions involving the ladder operators can be trans- 
lated into corresponding recurrence relation for the Laguerre polynomials and in particular we derive the Laguerre 
differential equation as follows. As the simplest example, consider the expression (PT|) . By using the relations among 

the radial eigenstates i?^'-* and namely, 

(X D\\l + 1,3-1) = r-'Dy\l + 1,3- 1) (36) 

it is possible to derive corresponding recurrence relation for the generalized Laguerre polynomials. By the change of 
variable = z and using the definition (l?71) . Eqn. ([55]) becomes 

^ lz-('+'^/2-i)e^Ae--z('+'^/2)i('+'^/2)(^). (37) 
3 dz 

Similarly, by using the normalized ladder operator Di, the recurrence relation 

\l + l,3-l)=Di\l,3), (38) 

in terms of Laguerre polynomials, becomes 

L^;+"^-'\z) = -^/;+t'^\z) (39) 
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where z = r . The recurrence relations (1371) and ((39)) can be used to obtain the second order differential equation 
that L^' ^^'^/^^ (^z) satisfy. By defining fi — I + d/2 — I for convenience and using the relations (|37)) and (p9)) one finds 

^^^.^(^) + (a^ + 1 - (^) + = (40) 



^(2.+i)W = -^e^ :r^-' H,,ir) (43) 



(^') = ^^(2.+i)M- (44) 



which justifies, a posteriori, that the introduced in ()27p corresponds to the generalized Laguerre polynomials. 

For d — 1 the quantum number corresponding to the angular momentum can take the value zero only while the 
radial quantum number becomes the principle quantum number for the 1-d harmonic oscillator and therefore for the 
values d = 1, / = the radial eigenfunctions become the energy eigenfunctions of the 1-d harmonic oscillator. By 
using the expression (|22l) for the radial eigenfunctions i?^'-* and the key functions i?^ ^ together with the definition of 
the Hermite polynomials, one has the relation 

i?2,(r) = (-l)^2^jU7i/2(,2) (41) 

where Hj(z) is the Hermite polynomial of order j and L'^iz) is a polynomial of order j in the variable z respectively 
[l6l |. Using this identification, it is easy to check the validity of the recursion relation (|37|) . This can be done by 
putting d ^ 1 and Z = in ([57)) . Thus, by the change of the variable r"^ — z, one finds 

i7^/^(z) = -z-VvAe-.V2iV2^^(,). (42) 

Using the relation ()4ip . and the recursion relation between the Hermite polynomials, namely, 

1 p d 
ypf dr 

which already follows from the construction above, one finds that 

Returning back to the recurrence relation (I37p . it can be can be put into a more useful form 

z-^e-^L^{z) ^ lAe-^z(^+i)L[;+i;(z) (45) 

with /i = / + d/2 — 1 for convenience. This form of the recurrence relation shows that ^ is the ladder operator between 
the functions e^'^ z^L^{z) and can be iterated k times to obtain 

L^{z) . ^^-e^|;e-^z('^+'=)i|;_\^}(z). (46) 

In terms of the radial coordinate and in terms of the radial eigenfunctions, this recurrence relation corresponds to the 
coordinate representation of the expression 

For fc = J, this yields the Rodrigues formula for the Laguerre polynomials rather then a recurrence relation among 
them. Similarly, the conjugate of the relation ()47p . namely, 

\l + k,j-k) D(i+k-i)D(i+k^2) ■ ■ ■ D(i+2)D(i+i)Di\l,j) (48) 

can be translated into the recurrence relation for the generalized Laguerre polynomials as 

i!rS(-)=(-i)'-|ii,"(-)- (49) 

The expressions for the radial eigenfunctions involving Z-changing operators and j) can also be translated 

into a recurrence relation for the generalized Laguerre polynomials. For example, using normalized ladder operators, 
the relation (I55| can used to derive recurrence relations for the generalized Laguerre polynomials as 

Lf-'\z)^z-^+^^^z^L';{z), (50) 

where ^ = I + d/2 — 1 and z = as before. Evidently, it is possible to derive further recurrence relations for the 
Laguerre polynomials from the expressions involving the ladder operators acting on the radial eigenfunctions. As the 
above typical examples illustrate, the normalization and recurrence properties of the Laguerre polynomials are built 
into the factorization method by construction. 
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IV. MORSE OSCILLATOR 

The above choice and the special way of writing the ladder operators -D(;+j) and for the harmonic oscilla- 

tor hamiltonian make it possible to obtain the useful properties of the generalized Laguerre polynomials. We will 
illustrate that the above approach also provides a framework for expressing energy eigenstates of factorizable (or 
supersymmetric) hamiltonians in terms of related special function. 

In this section we shall obtain the bound energy eigenstates of the Morse oscillator following the same approach as 
the harmonic oscillator. The Morse potential, which is important in the study of molecular vibrations, is given by 

V{x)^Voie~^^'= -26-^^^) (51) 

with the parameters Vq > and /3 having dimension inverse length with — oo < a; < oo [1, [l^liOl. It is also possible 
to transform the Morse potential problem into the two dimensional harmonic oscillator by an appropriate change of 



variables [21| . The finite number energy eigenstates for the bound states of the Morse potential can be expressed in 



units of length, momentum and energy of the Morse hamiltonian can be taken as hf3 and respectively. By 

defining the parameter A = tm\!i and shifting the coordinate by — -^InA, the Schrodinger equation for the Morse 
potential can be written as 



terms of the generalized Laguerre polynonmials by constructing appropriate ladder operators as follows. The natural 

md energy of the Morse hamiltonian can bi 
l^jT , and shifting the coordinate by — -g ' 

p2 + (e-- - A)'] ij ^ {X^ ~ \£\} ^ eV (52) 

which is convenient for our purposes. Therefore, the suitable ladder operators that factorize the Morse hamiltonian 
are of the form 

Dn = +ip-(e-"-n), (53) 
Di = _zp_(e--n). (54) 

As in the case of the harmonic oscillator, by using the same symbols for the operators, the coordinate representation 
of the conjugate ladder operators can be written conveniently as 

ax 

Di = -e"-+^"^^e-»---^, (56) 

by multiplying from the right with suitable integrating factors for the multiplicative terms and inverses of the in- 
tegrating factors from the left of the derivative. The recurrence relation for the ladder operators can be found to 
be 

=-01;A, +2n + l, n = 1,2,3... (57) 

by using the canonical commutation relation [x,p\ = i. It is convenient to work with the hamiltonian /i„ expressed in 
terms of the ladder operators (15^ and ([5^ , is given by 

hn =p^+ [e-^ - (n + i)]' = DiDn + {n+l). (58) 

The single key eigenket, which corresponds to the ground state, will be denoted by |n). The key eigenket is annihilated 
by the lowering operator. Thus D„\n) — yields the normalized eigenket 

= [r(2n)]V2 (59) 

corresponding to the ground state energy = n + j. Using the recurrence relation (|57|) and the definition of the 
Morse Hamiltonian ([55]) . one finds 

/^(«+.)^U) • • • ^1+1) = ■ ■ ■ ^1+1) l'^" + i2^'2{n + fc) + 1] I (60) 
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which determines all the excited states and the corresponding eigenvalues. In turn, by using the key functions ([59)1 
and the expression (|60p and the position representation of the raising ladder operator (|56p , the excited states can be 
written as 

with the corresponding energy s^n^j-^ = + X]fc=i[2(^i + fc) + 1]. The expression ([6T1) suggests that it is convenient 
to change to the variable z = 2e~^. Therefore, all the excited energy eigenstates can be written as 

in the new variable z. Note that with a suitable change of variables, the expression for excited states in terms of 
ladder operators (j62p resembles to the corresponding expression for the linear oscillator, where the excited states are 
obtained by repeated application of a single raising ladder operator on the ground state. Using the identity 

|-/W = (-lFp(^-+^^|-p<-^V(i) (63) 

where z = \/ p and for any continuous function /(z), it is possible to rewrite the expression (1621) in terms of the 
generalized Laguerre polynomials as 

^("-■)(^)-(f(ST7T^ (64) 

where z = 2e~^ and normalization constants of the eigenfunctions '0(i+j)(a;) can be found from (j35p . The eigenfunc- 
tions can easily be expressed back in terms of the parameters A and /3 of the Morse hamiltonian. In contrast to the 
d-dimensional harmonic oscillator, for the Morse hamiltonian there is only one quantum number yet the eigenfunctions 
are given in terms of the generalized Laguerre polynomials defined with two parameters. 



V. CONCLUSION 



The Schrodinger equations for isotropic harmonic oscillator and the Morse potential are studied using the factoriza- 
tion method of Infeld-HuU and for both of the potentials, the associated special function is identified as the generalized 
Laguerre polynomials and they are expressed in terms of appropriate independent variables which directly obtained 
from the factorization. In the factorization method of Infeld-Hull, the connection between the set of eigenfunctions 
and the special functions involved, is usually established, or rather verified, by comparing the results obtained from 
the algebraic methods with the results obtained by other methods \$\. In this respect, it is worth to emphasize that 
in the two cases presented above, the factorization is carried out without transforming the Schrodinger equation into 
a standard Hypergeometric-type differential equation by the change of the dependent and the independent variable 
at the outset and then employing the properties of the generalized Laguerre polynomials. In contrast, the generalized 
Laguerre polynomials are obtained directly using the expressions involving the ladder operators of the standard factor- 
ization method. Thus, we showed in the particular case of the harmonic oscillator that the feature of the Infeld-Hull 
factorization method which directly relates the energy eigenstates to the relevant special functions, that is, to the 
generalized Laguerre functions as in (|26p and also showed that any expression involving the ladder operators defined 
for the radial eigenfunctions in Section II can be used to define a corresponding recurrence relation for the generalized 
Laguerre polynomials with the help of the relation (|27p . 

Finally, we remind that the algebraic approach presented for the generalized Laguerre polynomials above is already 
known to apply to other special functions. A well-known example is the quantum mechanical angular momentum. 
Two different recurrence formulae for the associated Legendre polynomials can be obtained from the properties 

of the ladder operators Lip using the key functions P;^'(0)e^*''^ and the zero angular momentum states. By writing 
the coordinate representation of the ladder operators in a form similar to those of the harmonic oscillator, it is 
possible to find the Rodriguez type formulae from both oc \l,m) and (L+)"|Z,0) cx |;, m) for the P,"(6')'s. 

Another well-known example is the treatment of the quantum mechanical free particle in spherical coordinates. By 
omitting the term in the harmonic oscillator hamiltonian in three dimensions, the factorization method can be 
used to obtain the spherical Bessel and Neumann functions as well as their recurrence formulae. The different feature 
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of the latter example is that the normalization of spherical Bessel and Neumann functions (which involve (5-function 
normalization) can be shown to follow from the properties of the key functions corresponding to the zero angular 
momentum state. 
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